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Abstract 


Land subsidence is one of the most important geomorphic hazards that has a slow motion but 
destructive effects in the long run. For the spatial analysis of the amount of displacement, the 
data on changes of the ground surface from Sentinel Al for a period of 3 years (2017-2020) 
through the radar interferometric method and SNAP software were used. Field observations 
were also made to evaluate the research results. The results showed that the rate of land 
subsidence has increased from 2017 to 2020, from 2.6 cm in 2017 to 7.8 cm in 2020. Moreover, 
the groundwater hydrograph in the central plain of Ghaen has significantly decreased, so that 
the rate of water level drop in the plain level has been about 1.30 meters during 5 years (from 
2014 to March 2020), ie, about 26 cm per year. In addition, the analysis of spatial 
autocorrelation and Moran index (0.984) confirmed the clustering of subsidence event under the 
influence of water level changes in the study area. The correlation results also showed that there 
was a small correlation between the factor “changes in water level” and “land subsidence rate" 
(p = -0.138), which is an indirect relationship, i.e., the higher the groundwater abstraction (or 
the lower the water level), the higher the amount of subsidence. But due to the level of 
significance (sig = 0.585) this relationship is not significant. 
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حغرافیا و مخاطرات محیطی. سال یازدهم. شمارۀ چهل و eem‏ زمستان ۱ صص ۹4-٥‏ 
مقاله پژوهشی 


ویژه نامه (جالش جهانی فرونشست زمین: مدیریت بحران يا بحران مدیریت) 


تحلیل فضایی فرونشست سطح زمین با استفاده از تداخل سنجی راداری 


(موردمطالعه: دشت مر کزی شهرستان قاین) 
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حکیده 


فرونشست زمین یکی از مخاطرات مهم ژئومورفیک است که دارای حرکتی XS‏ بوده و در بلندمدت آثار 
مخربی برجای می گذارد. در این تحقیق به‌منظور تحلیل فضایی میزان جابجایی از داده‌های سنتینل 1 در 
بازه زمانی ۳ JUS‏ (۲۰۲۰- ۲۰۱۷) و روش تداخل سنجی راداری و نرم‌افزار (SNAP)‏ تغییرات سطح 


زمین را در بازه زمانی مذکور استخراج کرده‌ايم. همچنین مشاهدات میدانی برای صحت سنجی و بررسی 


نتایج تحقیق صورت گرفت. نتایج OLS‏ داد که میزان فرونشست زمین از سال ۲۰۱۷ تا سال ۲۰۲۰ افزایش 
داشته است؛ به‌طوری که از SA‏ سانتی‌متر در سال ۷ به VA‏ سانتی‌متر در سال YoYo‏ رسیده Soll‏ 
به‌طوری که Ol y‏ افت تراز آب در سطح دشت حدود ۱,۳۰ te‏ در طول ۵ سال (از سال ۲۰۱۶ تا مارس 


(Te‏ بوده است؛ cM‏ حدوداً سالی ۳۹ سانتی متر افت سطح آب‌های زیرزمینی را داشته‌ايم. ضمن اینکه 
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۱.۰ جغرافیا و مخاطرات محیطی شمارة چهارم 


تحلیل خودهمبستگی فضایی و شاخص موران (C AE)‏ خوشه‌ای بودن رخداد فرونشست تحت تأثیر 
تغییرات سطح آب در منطقه موردمطالعه را تأیید کرد. نتایج همبستگی نیز نشان داد که بین عامل تغییرات 
تراز سطح آب و نرخ yj‏ ت زمین» همبستگی کمی وجود داشته (0.138- (P=‏ که این ارتباط 
غیرمستقیم است؛ یعنی هر چه میزان برداشت آب زیرزمینی بیشتر باشد (یا سطح تراز آب کمتر شود)» 
میزان فرونشسست نیز بیشتر می‌شود؛ اما با توجه به سطح معناداری )0.585 (sigo‏ این ارتباط معنادار 


کلید واژه‌ها: ژئومورفیک. سطح تراز آب زیرزمینی» تداخحل سنجی» SNAP‏ قاین. 


Andie -۱ 


فروندٌ ت به‌عنوان یکی از وقایع ژئومورفیک: شامل نشست تدریجی و b‏ ناگهانی سطح زمین است که تحت 
تأثیر عوامل طبیعی و همچنین به‌واسطه دخالت‌های انسانی روی می‌دهد. 

اثرات فرونشست به شکل شکاف‌هایی بر روی سطح زمین» فرو چاله‌هه نشست منطقه‌ای و تشدید سیلاب 
شدگی نمایان می‌شود (مریخ پور و همکاران. .)۱۳٩۱‏ فروند ت طبیعی به‌واسطه فعالیت‌های مختلفی نظیر 


معدنکاری, برداشت آب‌های زیرزمینی و یا استخراج نفت» روندی تشدید شده به خود می‌گیرد تا جایی که به یکی از 


مخاطرات مهم ژئومورفولوژی تبدیل می‌شود. ایران به‌عنوان یکی از مهم‌ترین کشورهای مواجه با فرونشست تشدید 
asl‏ در آینده‌ای نه‌جندان دور با خسارات مستقيم 3 غیرستفیم جبران‌ناپذیری مواجه > Jal‏ شد. دراین‌بین» مناطق 
می‌باشند (جترسیماب و iU‏ ۱ گرچه در کوتاه‌مدت ممکن است در مناطق در حال فرونشست. خرابی به 


میزان گسترده مشاهده نود و حتی آثار سطحی حاصل از آن نیز به‌راحتی قابل تشخیص نباشد. اما بااین‌وجود در 


هستند؛ ازاین‌رو برای تشخیص میزان 9 9 23 ت» اندازه‌گیری و پایش تغییر شکل سطح زمین از اهمیت ویژه‌ای 


فرونشینی و نرخ فرونشست در LW‏ مختلف یک منطقه؛ (Y)‏ شناسایی عوامل وقوع فرونشست؛ (۳) ارائه راهکارهای 
مناسب به‌منظور کاهش نرخ فرونشست و سپس مدیریت و جلوگیری از خسارات ناشی از آن (افضلی و همکاران؛ 
۲ اسدزاده و همکاران ۱۳۹۵). 

برای اندازه‌گیری فروند ت زمین از روش‌های مختلفی ازجمله استفاده از داده‌های (GPS‏ ترازیابی دقیق و 


تکنیک‌های سنجش ازدور ب خصو ص تداحل سنچی راداری Y‏ استفاده می‌شود. des‏ سنچی راداری با 53 dow‏ 


1 GPS (Global Positioning System) 
2 Radar Interferometry 


سال یازدهم تحلیل فضایی فرونشست سطح زمین با استفاده از .... ۱۱ 


مصنوعی ' یکی از روش‌هایی است که در دو dao‏ اخیر انقلاب بزرگی در زمینه مطالعات تغییر شسکل پوستة زمین 
ایجاد کرده است (قابریل » AS‏ گر چه تداعل سنجی راداری اغلب به‌عنوان یک تکنولوژی جدید برای پایش 
دگرریختی‌های سطح زمین شناخته شده لیکن پتانسیل کامل کارایی‌های آن به‌وسیلة مجامع علمی به‌طور دقیق 
مشخص نشده و نتایج حاصل ازاین‌روش متعدد و پیچیده است (فررتی" و همکارانء ۲۰۰۷). تکنیک تداخل سنجی 
راداری با داشتن ویژگی‌هایی همچون پوشش زمینی پیوسته و قدرت SSE‏ مکانی و زمانی مناسب به‌عنوان یکی از 
بهترین روش‌های سنجش زدور در بررسی تغیبرات سطحی زمین شناخته شده است (روت و رقنل 4 Y V‏ در 
تداخل سنجی راداری» فاز تصاویری که در زمان‌های مختلف از یک منطقه برداشت شده‌اند. پیکسل به پیکسل مقایسه 
و سپس از هم کم می‌شود (فررتی و همکاران, ۲۰۰۷). 

در حال حاضر این روش در Oke‏ سایر روش‌های اندازه‌گیری مستقیم و غیرمستقیم. به‌عنوان کارآمدترین روش 
برای اندازه‌گیری تغیبرات سطح زمین با دقت و قدرت تفکیک مکانی بسیار UL‏ به شمار می‌رود (محمدی منش ۵ و 


همکاران. ۲۰۱۸). همچنین تا حد زیادی پی‌جویی فروند ت زمین در مناطق شهری را تسهیل کرده است 
V ULIS)‏ و همکاران. THVT‏ ازجمله مزایای دیگر این روش می‌توان به پوشش وسیع» نداشتن نیاز به کار میدانی؛ 
مقرون‌به‌صرفه بودن و امکان دستیابی به اطلاعات در هر شرایط آب و هوایی و دسترسی به آن با دقتی معادل با دقت 
اندازهگیری‌های سیستم مکان‌یابی جهانی و ترازیابی دقیق اشاره کرد. اولین تلاش‌ها برای نقشه‌برداری از فرونشست 
زمین در مساحت‌های بسیار بزرگ و در دوره‌های زمانی طولانی در جین صورت گرفته است (ژانگ" و همکاران, 
od MR‏ ماه کر تفر مه یکت انس ی کے PAL 13۱5۹۲ laa.‏ دیش 
قرار گرفته سه مجموعه داده ISAR‏ چهار ماهواره شامل تصاویر SAR‏ 86-1/2 از ۱۹۹۲ تا ۲۰۰۰ تصاویر 


SAR‏ از ۲۰۰۳ تا ۲۰۱۰ و تصاویر Y-RADARSAT‏ از سال ۲۰۱۲ تا ۲۰۱۶ مورداستفاده قرار گرفته است. 
درنهایت با بیش از ۱۲۰ اندازه‌گیری» سه نقشه از میزان فرونشست مربوط به سه دوره در منطقه مذکور AJ gi‏ شد که 
میزان فرونشست AN‏ میلی‌متر در سال (طی دوره ۱۹۹۲-۲۰۰۰ EM‏ میلی‌متر در سال (۲۰۰۲-۲۰۱۰) و ۵,۶ میلی‌متر 
در سال (۲۰۱2-۲۰۱۲) به دست آمد. فرونشست زمین ناشی از استخراج آب‌های زیرزمینی در منطقه شهری پکن با 
استفاده از تداخل سنجی راداری oF‏ مقدار فرونشست تا بیش از ۱۰۰ میلی‌متر در سال را در مناطق شرقی این شهر 
تأیید می‌نماید (چن" و همکاران (Y VV‏ در این مطالعه از Small Baseline NSAR S‏ برای پردازش تصاویر 


1 SAR Interferometry 
2 Gabriel 

3 Ferretti 

4 Rott, and Nagler 

5 Mohammadimanesh 
6 Castellazzi 

7 Zhang 

8 Chen 


Y‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


Envisat ASAR‏ به دست آمده بین سال‌های ۲۰۰۳ تا ۲۰۱۰ و تصاویر ماهواره‌ای TerraSAR-X‏ جمع آوری شده از 
سال ۰ نا ۱ برای بررسی فرونشست زمین استفاده شده است. علاوه بر این بین نقشه‌های نرخ فرونشست 
اینسار که از دو مجموعه osla‏ مختلف Envisat Jeol pea)‏ و (TerraSAR-X‏ حاصل می‌شود. همبستگی بالایی 
مشاهده شد. 

شارما و همکاران (VO‏ در مطالعه‌ای بر جزیره شرمن واقع در دلتای seal Slo‏ = سان واکین که یک منطقه pe‏ 
شهری با تغییرات کم زمانی است. نتایج حاصل از کاربرد تکنیک تداخل سنجی راداری برای اندازه‌گیری فرونشست 
آن منطقه را با استفاده از اطلاعات باند LL‏ طول‌موج ۲۳,۰۸ سانتی‌متر ماهواره UAVSAR‏ انجام دادند. در این مطالعه 


از داده‌های جمع‌آوری‌شده از ژانویه ۹ تا اوت ۶ برای ارزیابی تغییرات زمانی و مکانی فرو جزیره 
شرمن استفاده شد و نرخ فرونشست به‌صورت مقادیر متغیر بین ۵-۰ سانتی‌متر در سراسر محدوده به دست آورده 


شد. همچنین میانگین نرخ فرونشست منطقه تا حدود m‏ ۱,۳ سانتیمتر بر آوود شد. نتایج این مطالعه نشان داد که 


چندین عامل ازجمله نوع خاک عمق آب زیرزمینی» نوع کاربری زمین‌ها و ارتفاع زمین بر فرونشسست منطقه مؤثر 
است (شارما" و همکاران (Y V‏ مطالعات در نواحی مختلف ایران نشان می‌دهد در بیشتر دشت‌ها پدیده فرونشست 
در چند دهه اخیر فزونی یافته است. نرخ فرونشست در دشت اردبیل حدود ۸۰ تا ۱۵۶ میلی‌متر در سال محاسبه شده 
است (خلیفی و همکاران» (YAT‏ در مطالعه دیگری با استفاده از تکنیک تداخل سنجی راداری» بیشینه دامنه 
فرونشست دشت قزوین حدود ۳۰ تا ۳۵ میلی‌متر در سال ذکر گردیده است (بابایی و همکاران ۱۳۹۵). 

احمدی و همکاران (۱۳۹۷) در پژوهشی به مطالعه فرونشست دشت خرمده با استفاده از تکنیک تداخل سنجی 
راداری و بررسی مخاطرات of‏ پرداختند. برای این منظور از تصاویر ماهواره ENVISAT‏ اخذ شده در بازه زمانی 
۰۰۵-۳ و تصاویر \Sentinel S‏ اخذ شده در بازه زمانی ۲۰۱۶-۲۰۱۷ استفاده شد. میدان جابجایی» بیانگر پیشینه 
تغییر شکل برابر با ۳۵ میلی‌متر بر سال» در راستای دور شدن از ماهواره را نشان می‌دهد. نتایج به دست آمده از 


A‏ = با استفاده از سیستم اطلاعات مکانی (GIS)‏ با لایه‌های کاربری اراضی نیز جاه‌های بهره‌برداری در منطقه 


موردمطالعه نشان دهنده بالا بودن تراکم کاربری اراضی کشاورزی و نیز چاه‌های بهره‌برداری در منطقه است. همچنین 
با استفاده از تحلیل‌های مکانی. مخاطرات احتمالی فرونشست در منطقه موردمطالعه نیز بیانگر تأثیرپذیری قسمتی از 
مسیر خطوط راه‌آهن, راه‌های اصلی و بزرگراه از فرونشست است. 

تورانی و همکاران (۱۳۹۷) با Alas‏ فرونشست در شهر گرگان با استفاده از روش تداخل سنجی راداری دریافتند 
که در دوره زمانی بین بهمن ۱۳۸۵ تا آذر ۱۳۸۸ محدودۀ شمال گرگان متحمل فرونشست شده است. این محدوده 


تقریباً روند شرقی- غربی دارد و میزان فرونشست در Ol‏ چیزی حدود ۸ سانتیمتر محاسبه شد. نمودارهای تراز 


1 Sharma 


سال یازدهم تحلیل فضایی فرونشست سطح زمین با استفاده از .... yey‏ 


سطح آب و میزان بارندگی گرگان در سال‌های ۱۳۸۵ تا ۱۳۸۹ علیرغم نوسانات فصلی روند نزولی OUS‏ می‌دهند. 
پیامد این افت سطح آب زیرزمینی می‌تواند منجر به فرونشست زمین در شهر گرگان شده باشد. 

مهرابی (۱۳۹۷) در پژوهشی به شناسایی شواهدی بر وجود گنبد نمکی مدفون و جدید در ناحیه زاگرس با 
استفاده از روش تداخل‌سنجی تصاویر راداری ستتینل-۱ و ایسار پرداخت. بدین منظور از تصاویر راداری ایسار 
(ASAR)‏ مربوط به سال‌های ۲۰۰۷ و ۲۰۱۲ و ستتینل-۱-۹60/061(۱) مربوط به سال‌های ۲۰۱۶ و ۲۰۱۷ استفاده 
شد. جهت تعیین میزان تغییرات پوسته زمین ضمن انجام فرایند واپیچش فاز بر روی تصاویر, تداخل نگاشت هر یک 
از تصاویر استخراج شد. با توجه به نتایج به دست آمده از تحلیل تصاویر ایسار نرخ رشد گنبد نمکی مذکور بین 
سال‌های ۲۰۰۷ تا ۲۰۱۲ به میزان ۱,1 سانتی‌متر در سال بوده است. طبق تصاویر تداخل نگاشت به دست آمده از 
تصاویر سنتینل-۱ سرعت صعود این گنبد نمکی در بین سال‌های ۲۰۱۶ تا ۲۰۱۷ افزایش داشته است و به ۲,۹ 
سانتی‌متر در سال رسیده است. نتایج تحقیق نشان داد که گنبد نمکی احتمالی فعال است. 

مقصودی و همکاران (۱۳۹۸) در پژوهشی به بررسی رفتار فرونشست زمین در منطقه غرب تهران با استفاده از 
تصاویر سنجنده سنتینل-۱ و تکنیک تداخل سنجی راداری مبتنی بر پراکنش گرهای دائمی پرداختند. در این مقاله یک 
سری زمانی دو ساله شامل ۳۰ تصویر سنجنده ستتینل- ۱ با استفاده از روش تداخل سنجی راداری مبتنی بر 
پراکنشگرهای دائمی در دو ناحیه مطالعاتی مورد پردازش قرار گرفت. نتایج» جابجایی سالانه ۱۱ سانتی‌متری در 
جنوب منطقه موردمطالعه را نشان cols‏ اما جابجایی تنها محدود به مناطق دشتی و غير شهری نبوده بلکه به مناطق 
صنعتی و شهری نیز گسترش HL‏ است. ارزیابی نسبی نتایج در منطقه اول با استفاده از داده‌های مسیر پایین گذر ۳۱ 
انجام گرفت که با توجه به احتلاف ناچیز برآورد دو مسیر ۷۹ و ۳۱ نواحی مذکور جابجایی افقی بسیار کمی داشته و 
oles‏ بردار جابجایی را حرکت قائم تشکیل می‌دهد. همچنین. جابجایی قائم با استفاده از دو مدار بالاگذر و پایین گذر 
برای منطقه مطالعاتی دوم برآورد گردید که با بهره‌گیری از داده‌های یک ایستگاه GPS‏ صحت سنجی شد. 

_آروین و همکاران (YA)‏ با استفاده از تصاویر ماهواره‌ای راداری ۱-۹6۳6 سال‌های ۱۳۹۳ و ۱۳۹۷ و 
روش تداخل‌سنجی تفاضلی راداری و سیستم اطلاعات جغرافیایی به مدل‌سازی مکانی فرونشست زمین در جنوب 
حوزه آبخیز میناب پرداختند. qld‏ پردازش‌های انجام‌شده بر روی این تصاویر نشان داد که در دوره مطالعاتی» حدود 


۱۳ سانتی‌متر فرونشست داشته است. بعد از صحت سنچی shig‏ انجام‌شده با واقعیت‌های زمینی (شهریور «(AV‏ 


به‌منظور تشسخیص رابطه فرو با میزان تغییرات و افت place‏ آب زیرزمینی؛ ارتفاع سطح زمین و شیب 
cle ules‏ فضایی انجام و میزان همبستگی هر یک از فاکتورهای مذکور با تراکم رخداد فرونشست محاسبه گردید. 
به نظر می‌رسد از مجموع ۰ دشت کشون بیش از نیمی از آن‌ها در معرض تفست باشند. اگر زمانی این پدیده 


تنها در برخی استان‌ها مانند کرمان و یزد به وقوع می‌پیوست امروزه این پدیده در دشت‌ها و شهرهایی نظیر epe‏ 


Pt‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


تهران» کاشمر «QUSS‏ خراسان جنوبی و .. روی داده oul‏ (شریفی (OA MY. LS‏ این مطالعه سعی در تهیه نقشه 
که برای انجام آن, از داده‌های تصاویر ماهواره‌ای ستتینل ۸۱ در نرمافزار (SNAP)‏ استفاده شده است. 


۲- منطقه موردمطالعه 

دشت مرکزی قائنات در یک حوضه آبریز بسته بین ۵٩‏ درجه و ۳ دقیقه تا ٥۰‏ درجه و ۳ دقیقه طول شرقی و ۳۳ 
درجه و ۱ anis‏ تا ۳۳ درجه و ۵۶ دقیقه عرض شمالی واقع شده است (شکل C‏ دشت قائنات یک دشت بسته 
آبرفتی بوده که توسط رسوبات دوران چهارم به‌صورت کامل پوشیده شده است. ازنظر توپوگرافی سطح دشت نسبتاً 
ناهموار است که با یک شیب بسیار کم VA)‏ درصد تا ۳,۳ درصد) به نواحی کم ارتفاع شرق دشت ختم می‌شود. 
حداکثر ارتفاع در کوهستان‌های منتهی به دشت ۲۵۵۲ متر در جنوب Of‏ و حداقل آن ۸٩۰‏ متر در شرق دشت است. 
حوضه آبریز دشت مرکزی قاين در وسعتی معادل ۸ کبلومترمربع دارد که از این مقدار حدود ۰ کبلومترمربع 
خراسان جنوبی به‌ویژه در شرق استان» کاهش سطح آب‌های زیرزمینی» به‌عنوان منبع اصلی تأمین آب شرب شهر 
قاین و روستاهای محدوده موردمطالعه را به دنبال داشته است. برداشت بی‌رویه آب زیرزمینی باعث افت ممتد سطح 


آب و درنتیجه موجب ایجاد و گسترش فرونشست در استان و به‌ویژه در منطقه مطالعاتی شده است. 


59* 30'0"E 
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شکل ۱- موقعیت منطقه موردمطالعه در تقسیمات سیاسی کشوری 


سال یازدهم تحلیل فضایی فرونشست سطح Om)‏ با استفاده 5l‏ — ۱۵ 


۳- مواد و روش 
Y‏ -3 - داده‌های مورداستفاده در تحقیق 


روش کار در این تحقیق تحلیلی و میدانی است. در این تحقیق از داده‌های راداری که شامل تصاویر سنتینل ۱ 


جهت محاسبه نرخ فرونشست در بازه زمانی یک ماه برای سال‌های ۲۰۱۷ تا ۲۰۲۰ استفاده شده است و پس از 


محاسبه نرخ جابجایی سطح زمین برای منطقه موردمطالعه به‌منظور بررسی بیشتر SE‏ افت تراز آب زیرزمینی بر 


سعی شد که علاوه بر نوسانات تراز GT‏ پیزومترها؛ پراکندگی چاه‌های بهره‌برداری» میزان برداشت آب» 


فروزه 
نقشه پهنه‌بندی میزان تخلیه سالانه آب زیرزمینی نیز موردبررسی قرار بگیرد مطالعات تداخل سنجی راداری در این 
تحقیق توسط VE‏ تصویر سنتینل و به کمک نرم‌افزار SNAP‏ انجام گرفته است. برای کاهش نویز و افزایش همبستگی 
زمانی اینترفروگرام‌ها تا حد امکان سعی شده از زوج تصاویر راداری با بازه زمانی کوتاه (حدود ۰ روز) استفاده 
شود. ترکیب داده‌های راداری از مدارهای صعودی و نزولی نیز می‌تواند جهت بهبود مدل‌های رقومی زمین با به 
دست آوردن جابجایی‌ها در جهات مختلف و با استفاده از اینترفرومتری مورداستفاده قرار گیرد. با عنایت به اينکه در 
این تحقیق هدف پایش میزان جابجایی‌های سطح زمین است. تنها از داده‌های صعودی استفاده شده است؛ بنابراین 
اساس از تعداد VE‏ تصویر مربوط به دوره زمانی ۲۰۱۷ تا ۲۰۲۰ مورداستفاده قرار گرفت که در جدول Y‏ اطلاعات 


ثبت شده است. 


حدول ۱- مشخصات عمومی تصاویر استفاده شده در تداخل سنجی راداری 


sies path | Ascending/Descending | Polarization pc Frequency Mission‏ ین 
IW 86 Ascending VV-HV 16083 C-Band Sentinel-1‏ 2017/04/10 
IW 86 Ascending VV-HV 16958 C-Band Sentinel-1‏ 2017/06/09 
IW 86 Ascending VV-HV 17308 C-Band Sentinel-1‏ 2017/07/03 
IW 86 Ascending VV-HV 18533 C-Band Sentinel-1‏ 2017/09/25 
IW 86 Ascending VV-HV 18708 C-Band Sentinel-1‏ 2017/10/7 
IW 86 Ascending VV-HV 19933 C-Band Sentinel-1‏ 2017/12/30 
IW 86 Ascending VV-HV 20108 C-Band Sentinel-1‏ 2018/01/11 
IW 86 Ascending VV-HV 21158 C-Band Sentinel-1‏ 2018/03/24 
IW 86 Ascending VV-HV 21333 C-Band Sentinel-1‏ 2018/04/05 
IW 86 Ascending VV-HV 22558 C-Band Sentinel-1‏ 2018/06/28 
IW 86 Ascending VV-HV 22733 C-Band Sentinel-1‏ 2018/07/10 
IW 86 Ascending VV-HV 23783 C-Band Sentinel-1‏ 2018/09/20 
IW 86 Ascending VV-HV 23958 C-Band Sentinel-1‏ 201810/02 
IW 86 Ascending VV-HV 25183 C-Band Sentinel-1‏ 2018/12/25 
IW 86 Ascending VV-HV 25358 C-Band Sentinel-1‏ 2019/01/06 
IW 86 Ascending VV-HV 26583 C-Band Sentinel-1‏ 2019/03/31 
IW 86 Ascending VV-HV 26933 C-Band Sentinel-1‏ 2019/04/24 
IW 86 Ascending VV-HV 27808 C-Band Sentinel-1‏ 2019/06/23 
IW 86 Ascending VV-HV 27983 C-Band Sentinel-1‏ 2019/07/05 
IW 86 Ascending VV-HV 29208 C-Band Sentinel-1‏ 2019/09/27 


ER!‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


ei s path | Ascending/Descending | Polarization pes Frequency Mission‏ رل 
IW 86 Ascending VV-HV 29383 C-Band Sentinel-1‏ 2019/10/09 
IW 86 Ascending VV-HV 30433 C-Band Sentinel-1‏ 2019/12/20 
IW 86 Ascending VV-HV 30608 C-Band Sentinel-1‏ 2020/01/01 
IW 86 Ascending VV-HV 31133 C-Band Sentinel-1‏ 2020/01/06 


۲-۳- روش تحقیق 

برای پردازش تصاویر ماهواره‌ای در محیط نرم‌افزار 5۷۸ از روش تداخل سنجی راداری با دریچه مصنوعی و 
ee cs‏ مج adc Basle ul in SAR) cas‏ اسر jeu‏ سنج urls‏ قار تساو ال سلاو 
موقعیت‌های تصویربرداری و یا زمان‌های تصویربرداری مختلف پیکسل به پیکسل مقایسه گردید. از تفاضل‌گیری 
بین این مقادیر» تصویر اینترفروگرام تولید شد. 

پس از تشسکیل تداخل da lS‏ یک شبکه از تصاویر ایجاد شد و با استفاده از روش کمترین مربعات. مقدار 
جابجایی هر پیکسل تخمین زده شد. رابطه فاز ایتترفرومتری Geb‏ رابطه )1( محاسبه گردید. 

رابطه (۱) 

yi = ly, | exp(ig,) 


y: = lyzl exp(iq;) 
yia = lyilly;l exp(i(@: — 92) 


در این رابطه؛ Vi‏ و ya‏ سیکنال‌های دو پیکسل متناظر است که اولی تصویر اصلی (master)‏ و دومی تصویر فرعی 
Glave)‏ است. پس از ضرب مختلط این دو سیگنال رابطه دیگری به دست می‌آید که در این رابطه [yy Mya]‏ دامنه 
Pi — 92‏ و برابر فاز اینترفروگرام ias‏ در این رابطه تخییر فاز به‌اندازه TY‏ معادل جابجایی به‌اندازه صف طول‌موج 
مورداستفاده توسط ماهواره است که نشان‌دهنده یک فرینج کامل در اینترفروگرام است (اکبری و همکاران, ۱۳۹۹ 

فاز اینترفروگرام ایجاد شده ناشی از مولفه‌هایی نظیر ملفه مداری» توپوگرافی» جابجایی؛ اتمسفر و نویز است که 
هر یک از این پارامترها سیب تغییر فاز می‌شود. ارتباط اختلاف‌فاز ایجاد o s‏ در dels‏ سنجی راداری 3 موْلفه‌های 
ذکرشده به‌وسیله رابطه زیر تعیین می‌شود. 

رابطه (0) 


T A®peformation + At Atmosphere noise 


Ap, = O21 ¬ P2 = AtpGeometry T A®ropography 


سال یازدهم تحلیل oles‏ فرونشست سطح زمین با استفاده از .... 1۷ 


موقعیت نسبی تعدادی از عوارض زمینی در بازه زمانی بین دو تصویربرداری SAR‏ به‌صورت جزئی تغییر می‌کند. 
این تغییرات ناشی از پدیده‌هایی نظیر فرونشست. زمین‌لرزه» زمین‌لغزش و یا حرکت گسل‌ها است. این اثر تغییر فازی 
زا مسقل از خط ما ایجاد می کد که از رابطه زیر تین 3h ge‏ 

(Y) رابطه‎ 


4m 
Apa = 44d 


در این رابطه؛ d‏ تصویر جابجایی نسبی عارضه در راستای رنج مایل است. پس از مسطح سازی اینترفروگرام فاز 
تداخل سنجی شامل هر دو اثر ارتفاعی و جابجایی خواهد بود؛ بنابراین تداخل سنجی در این حالت مطابق رابطه زیر 
تعیین می‌شود (زندی و همکاران. ۱۳۹۹). 

(ads 


4n B 4m 
c اف‎ mes d 


۵ = A Rsin8 A 


تحلیل فضایی فرونشست سطح زمبن با استفاده از تداخل سنجی رادار ی 
سنطقه سورد سطالعه( دشت مر کزی قاننات» 


کی 


ED ER 


ph‏ راداری 
(نرم اخزار Snap‏ 


داده al>‏ جا AW‏ و الاح جس cca‏ فاز 


J (GIS-SPSS (نرم 3# ار‎ 


ahas XL A شتا‎ 
Nd فرونشست‎ — 


لس 


۱ تولہد نقشد نهایی فر وفشست 


شکل ۲- فلوچارت مراحل انجام تحقیق 


۱۸ جغرافیا و مخاطرات محبطی شمارة چهارم 
-٤‏ نتایج و بحث 

ME‏ پایش فرونشست زمین 

استخراج آب زیرزمینی و افت سطح آب‌های زیرزمینی به میزان فوق سبب کاهش فشار هیدرواستاتیک و خروج 
آب از خاک‌های سطح مخروطی دشت که مشتمل بر لایه‌های ریز و درشت است. منجر به تحکیم و کاهش حجم 
خاک و درنهایت رویداد نشست تدریجی زمین است؛ بنابراین علیرغم وجود علت‌های مختلف قابل‌مطالعه در ایجاد 
نشست؛ در تحلیل نمودار و نیز با تعیین افت سطح آب می‌توان صحت فرض این پژوهش را در دشت قائنات تأیید 
نمود. در این تحقیق با روش تداخل سنجی تفاضلی داده‌های راداری میزان و دامنه جابجایی سطح زمین برای دشت 
مرکزی CoU‏ تبیین شده است. برای منظور تولید نقشه‌های نشان دهنده الگوی فضایی میزان و دامنه نشست از زوج 
تصوير سنجنده‌های راداری در بازه زمانی ۲۰۱۷ تا ۲۰۲۰ (زمستان ۱۳۹۵ تا زمستان ۸ جهت انجام عملیات 


dels‏ سنجی استفاده شده است. 


بررسی الگوی فصلی فرونة ت زمین در این بازه زمانی مشخص کرد که میزان فرونشست زمین در دشت 
مرکزی قائنات. در فصول مختلف سال و پهنه‌های مختلف متفاوت بوده است. بر همین اساس در سال ۰۲۰۱۱۷ 
بیشترین میزان فرونشست (۰,۰۸۱۱ متر یا ۸ سانتی‌متر) مربوط به فصل بهار بوده که در نواحی شمالی دشت قائنات 
pe lel ont‏ کن مان قر ات مط هه فل وق SiO A La osa‏ ری دار 
ac‏ کو ead‏ مهار برای ۱۱ سای ی رھ Ned, oly yas‏ دنت قاقات dla‏ اقات ادت انا دز 


pL‏ فصول بیشترین فرونشینی در نواحی co gom‏ غربی دشت اتفاق افتاده است. در سال ۹ نیز ميزان فرونشست 


در حال افزایش است. در این سال نیز بیشترین ميزان فرونشست زمستان اتفاق افتاده که برابر با ۰,۱۹۹ متر (۱۹,۹ 
سانتی‌متر) بوده که در نواحی جنوب شرقی دشت قائنات اتفاق افتاده است. میزان فرونشست فصل تابستان نیز بالا و 


برابر با ۱۱,۵ سانتی‌متر بوده که در نواحی جنوبی دشت موردمطالعه اتفاق افتاده است. 


سال ehh‏ تحلیل فضایی فرونشست سطح ME‏ با استفاده از .... ۱۹ 


شکل ۳- الگوی فصلی فرونشست دشت مر کزی قائنات طی سال‌های ۲۰۱۷ تا ۲۰۲۰ 


مأخذ: یافته‌های پژوهش ۱۳۹۹ 


در ادامه روند سالانه فرونشست دشت قائنات طی دوره موردمطالعه نیز موردبررسی قرار گرفت. با توجه به شکل 
(E)‏ الگوی فرونشست برای منطقه موردمطالعه در جهت شمال غربی-جنوب شرقی قرار دارد. Gb‏ نقشه موردنظر 
بیشترین نرخ فرونشست در بازه زمانی یک‌ساله برای سال ۱۳۹۷ حاصل از داده‌های ۰,۰۰۲,۹۵00061 متر SY‏ 


11۰ جغرافیا و مخاطرات محیطی شمارة چهارم 


سانتیمتر) و در قسمت جنوب شرقی منطقه موردمطالعه است. همچنین میزان ۰,۰۸۳ متر (۰,۸۳ سانتیمتر) بالاآمدگی 
سطح زمین در بازه زمانی موردنظر در قسمت غرب و شمال غربی منطقه موردمطالعه نیز برآورد شده است. 

مقدار قرو تتت Ys Y E PAV E Lo‏ ور خال Ala‏ اس که از YS‏ در« مر در ال ۲۰۱۷ به HO‏ 
در سال ۲۰۱۸ ٠,۰۵٤۲‏ متر در سال ۲۰۱۹ و در سال ۲۰۲۰ به ۰,۰۷۸ متر VA)‏ سانتی‌متر در سال) می‌رسد. پهنه‌های 
در معرض فرونشست سالانه بیشتر منطبق بر نواحی جنوب شرقی دشت قائنات یعنی نحای پیرامونی شهر اسفدن 


است اما نواحی شمال غربی دشت و نواحی پیرامونی شهر قاین نیز با فرونشست همراه است. 
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-Y-i‏ شرایط هیدروژئولوژی منطقه 
برای مشاهده ارتباط فرونشست با برداشت آب زیرزمینی در منطقه هیدروگراف‌های معرف آب زیرزمینی دشت 


مرکزی شهرستان قائنات asl jl‏ شده است. 


سال یازدهم تحلیل فضایی فرونشست سطح زمین با استفاده از .... ۱۱۱ 


با توجه به وجود شواهدی مانند: ۱) بردارهای GPS‏ و مشاهده میزان تغییر شکلی که این بردارها نمایش می‌دهند. 
۲ حذف تصاویری نزدیک به تاریخ وقوع زمین‌لرزه که عامل اصلی جابه‌جایی در بازه زمانی خیلی کوتاه هستند و ۳) 


تطابق بیشترین نرخ فروند 7 با محدوده آبخوان‌های آبرفتی. می‌توان گفت که عامل اصلی رخ دادن فرونشست 


زمین در منطقه موردمطالعه» برداشت بی‌رویه منابع آب زیرزمینی و افت شدید تراز آب است. نشستی که تحت تأثیر 
پدیده تحکیم بر اثر خروج آب از ابه‌لای ذرات خاک‌های ریزدانه (رس و سیلت) حاصل می‌شود و وقوع Ol‏ تدریجی 
و تابع زمان است را نشست غير الاستیک گویند (رحیمیء ۱۳۸۶)؛ یعنی در این حالت بر خلاف نشست الاستیک؛ 
تغییر شکل بلافاصله بعد از اعمال نیرو به وقوع نمی‌پیوندد. بلکه تدریجاً و با پیشرفت زمان رخ می‌دهد. واحدهای 
رسی و سیلتی به علت ریزدانه بودن» تخلخل UL‏ و ضریب نفوذپذیری قابل‌ملاحظه نقش SAIS‏ در مناطق فرونشست 
بر عهده دارند؛ زیرا نشست حاصل از برداشت آب از این لایه‌ها غیرقابل بازگشت و دائمی خواهد بود (ون جی " و 
همکاران» ۲۰۱۳) و سبب کاهش ضریب ذخیره آبخوان و در موارد بحرانی سبب ابود شدن آبخوان برای هميشه 
می‌شود. چاه‌های بهره‌برداری در منطقه مهم‌ترین عامل تخلیه آب از دشت‌ها هستند و ازانجاکه هر چاه با تخلیه آب 
منفذی محدوده شعاعی اطراف خودش سبب کاهش فشار آب منفذی می‌شود تعادل نیروها به هم می‌خورد و تنش 
مؤثر در محدوده اطراف چاه افزایش می‌یابد. برای جبران کاهش فشار آب منفذی آبخوان متراکم می‌شود و این امر به 
نشست غیرالاستیک در لایه‌های ریزدانه می‌انجامد. به‌منظور بررسی ارتباط فرونشست با برداشت آب زیرزمینی در 
منطقه» در این بخش هیدروگراف‌های معرف آب زیرزمینی دشت مرکزی قائنات به همراه تغییرات تراز بلندمدت آب 
زیرزمینی پیزومترهای مختلف ارائه خواهد شد زیرا رفتار پیزومترها به‌خوبی بیانگر روند کاهش فشار آب منفذی در 
آبخوان است. اندازه‌گیری تراز آب. از مسطح زمین تا سطح آب موجود در چاه صورت می‌گیرد و بنابراین» هرچه 
میزان بیشتر باشد. در حقیقت افت سطح آب چاه را OUS‏ می‌دهد. در شکل ۵ مشاهده می‌شود که هیدروگراف معرف 
آب زیرزمینی در دشت مرکزی قائنات افت محسوسی در طی چندین سال اخیر داشته به‌طوری که عمق آب‌های 
زیرزمینی به میزان ۱,۳۰ متر افزایش یافته است. این مسئله بیانگر این است که میزان افت تراز آب در دشت حدود 
۳۰ متر در طول ۵ سال (از سال ۱۳۹۳ تا اسفند ۱۳۹۸) بوده است» یعنی حدوداً سالی YA‏ سانتی‌متر افت سطح 
آب‌های زیرزمینی را داشته‌ایم. 

هیدر وگراف واحد منطقه نیز که با استفاده از داده‌های مربوط به چاه‌های مشاهده‌ای طی دوره موردمطالعه ترسیم 
شده است. اگرچه چندین بالاآمدگی سطح تراز آب زیرزمینی را (ازجمله در آذر ۱۳۹۵ تا خرداد OLE (PAT‏ می‌دهد 
اما روند مسینوسی نمودار حکایت از کاهش قابل‌توجه سطح آب زیرزمینی است. همان‌طور هم که در شکل ۵ 


مشسخصی است روند افزایش سالانه برداشت آب تا ماه gai gd‏ و مهر ادامه دارد و از مهرماه پا OUI‏ که بارندگی‌ها 


1 Wen-jie 


NY‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


افزایش یافته» مجدد میزان برداشت آب با روند کاهشی روبه‌رو شده است که دلیل این امر هم می‌تواند کاهش She‏ 
آبی» شروع بارندگی‌های پاییز و زمستان و توقف روند برداشت آب از چاه باشد. بر اساس بررسی انجام شده طی 
دوره موردمطالعه سطح تراز آب زیرزمینی در حال افزايش است و سال ۱۳۹۸ بیشترین افزایش سطح تراز آب (افت 


سطح (OT‏ را داشته است همچنین بیشترین شیب مشاهده شده مربوط به اسفند WAV‏ تا فروردین ۱۳۹۸ است. 


1393 1394 1395 | 1396 | 1397 1398 


شکل ۵- میزان برداشت آب از چاه‌های نمونه واقع در دشت مرکزی قائنات از فروردین ۱۳۹۳ تا اسفند ۱۳۹۸ 


شکل V‏ روند افزايش عمق دسترسی به آب زیرزمینی برای هرماه ارائه شده با توجه به این شکل مشاهده می‌شود 
که بیش‌ترین مقدار روند کاهشی در ماه مهر و شهریور و کم‌ترین مقدار در ماه فروردین است. به‌طوری‌که از شکل نیز 
می‌توان استنباط کرد ale‏ شیب‌های خط روند تغییرات تراز آب زیرزمینی در تمام ماه‌های سال منفی و ما بین ۰,۲۱ تا 
۹ متر در سال تغییر می‌کند؛ به‌عبارت‌دیگر سطح ایستابی آبخوان دشت قائنات به‌طور متوسط در طول T‏ سال بین 
۱ در فروردین‌ماه و ۰/۳۹ متر در آبان کاهش پیدا می‌کند؛ که یکی از عوامل افت سطح آب زیرزمینی را بهره‌برداری 
غیراصولی از چاه‌های حفر شده می‌توان معرفی کرد. 


سال يازدهم تحلیل فضایی فرونشست سطح Om)‏ با استفاده از .... ۱۱۳ 


شکل *- روند تغییرات تراز آب زیرزمینی دشت مر کزی قائنات از سال ۲۰۱۶ تا سال ۲۰۲۰ 


مأخذ: یافته‌های پژوهش ۱۳۹۹ 


برای این که ملاحظه شود چه رابطه‌ای بین میزان برداشت آب و ارتفاع برقرار بوده ضریب همبستگی بین 
متغیرهای مذکور موردمحاسبه قرارگرفته؛ به‌این‌ترتیب که ضریب همبستگی بین ارتفاع به‌عنوان متغیر مستقل و سطح 
برداشت آب در سطح اشنا موردمطالعه به‌عنوان متغیر وایسته محاسبه‌شده AV SX) Cul‏ ضریب بەدستآمده 
only‏ اساس براير با ۶ بوده که همبستگی ضعیف و معنی‌داری را بین دو متغیر مذکور نشان می‌دهد؛ 


به‌عبارت‌دیگر در نواحی با ارتفاع بیشت افزایش سطح تراز آب زیر زمینی نیز به نسبت بیشتر بوده است. 


y —0.5816x + 1241.5 
R?= 0.0271 


Water harvesting rate 


)۱۳۹۹ میزان همبستگی بین متغیرها (مأخذ: یافته‌های پژوهش.‎ -V JRE 


Mí‏ جغرافیا و مخاطرات P‏ شمارة چهارم 


مشخصات آماری چاه‌های موردمطالعه و میزان میانگین سالانه تراز آب در دوره آماری سال‌های ۲۰۱۷ تا ۲۰۲۰ 
برای چاه‌های موردبررسی در جدول ۲ درج شده است. همان‌طور که مشاهده می‌شود ایستگاه یا چاه پهنایی با 


میانگین افت سطح تراز ۹ متر در طول دوره موردمطالعه» شدیدترین روند منفی را داشته است 4$ در غرب 
منطقه موردمطالعه واقع شده است. 


حدول ۲- متوسط سطح تراز آب زیرزمینی آبخوان آبرفتی جاه‌های نمو 43 در محدوده موردمطالعه 
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مأخذ: سازمان آب منطقه‌ای استان خراسان جنوبی 


سال یازدهم تحلیل فضایی فرونشست سطح زمین با استفاده از .... Mo‏ 


همچنین نتایج آزمون همگنی روند تراز آب زیرزمینی دشت موردمطالعه نشان می‌دهد که مقدار 62(محاسبه شده 
برای ایستگاه‌ها بزرگ‌تر از 262 متناظر جدول با درجه آزادی ۱۷ است؛ بنابراین» ایستگاه‌ها از نظر داشتن روند تغییرات 
تراز آب زیرزمینی در سطح ۵/ همگن نیست. نتایج به دست آمده در شکل زیر قابل‌مشاهده است؛ به‌طوری‌که چاه 
پهنایی» شرق روستای خطیبی و مسیر راه سیستانک بیشترین افت سطح تراز آب را داشته و چاه‌های شمال اسفشاد و 


مسیر راه گرماب کمترین میزان افت را در سطح تراز آب داشته است. 
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شمال اسفدن جنوب شرق اسفدن 


شکل -A‏ مقایسه میزان افت سطح تراز آب در چاه‌های موردمطالعه 


مأخذ: یافته‌های پژوهش» ۱۳۹۹ 


شکل ٩‏ ارائه شده است. 
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شکل 4- مشاهدات میدانی نویسندگان از بررسی فرونشست دشت مر کزی قائنات 


۳-۶- تأثیر هیدروژئولوژی منطقه بر روی فرونشست 


به‌منظور درک بهتر تأثیر هیدروژئولوژی منطقه بر روی فرون" ت» نتایج حاصل از تداخل سنجی راداری 
به‌صورت نمودارهای سری زمانی ترسیم شد و با نمودارهای تغییرات سطح آب زیرزمینی در پیزومترهایی که در 
همان نقاط قرار داشتند» مقایسه شد (S)‏ منحنی‌های آبی سطح آب زیرزمینی را بر حسب متر و منحنی‌های 
نارنجی نرخ میانگین فرونشست سالانه را بر حسب سانتی‌متر بر سال در کل منطقه موردمطالعه OUS‏ می‌دهند. بر طبق 
این مقایسه مشاهده می‌شود که روند فرونشست. تطابق قابلملاحظه‌ای با روند کاهش سطح آب زیرزمینی در مناطق 
مختلف داشته و در مناطقی که فرونشست بیشتری رخ داده است. بیشترین میزان افت سطح آب زیرزمینی نیز مشاهده 
می‌شود؛ اما نکته قابل توجه این بوده که بیشترین میزان فرونشست در مناطقی رخ داده است که علاوه بر افت شدید 
تراز آب. لایه‌های ضخیم رس و سیلت با قابلیت فشردگی دائم بالا نیز وجود دارند. همچنین در برخی مناطق با وجود 
افت شدید تراز آب» هیچ نشست قابل‌ملاحظه‌ای رخ نداده است که بررسی اطلاعات دانه‌بندی لایه‌های زیرین نشان 


سال یازدهم تحلیل فضایی فرونشست سطح زمین با استفاده از .... ۱۷ 


نکته قابل توجه این است که از سال ۷ به بعد تراز آب زیرزمینی افت شدیدی داشته که این امر با نرخ 
فرونشست بیشتر در منطقه همراه است. دلیل این امر به cla JL See‏ شدید سال‌های اخیر مربوط می‌شود که منجر 


به افزايش تقاضا برای آب و اضافه برداشت شدید آب زیرزمینی در منطقه شده است. 
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شکل ۰- روند سالانه و فصلی سطح آب زیرزمینی و فرونشست زمین در دشت قائنات (سال‌های ۲۱۷ تا 


(T° 


مأخذ: یافته‌های پژوهش ۱۳۹۹ 


تحلیل خودهمبستگی فضایی و شاخص موران نیز خوشه‌ای بودن رخداد فرونشست تحت تأثیر تغییرات سطح 
آب در منطقه موردمطالعه را تأیید کرد. خودهمبستگی به رابطه بین مقادیر باقیمانده در طول خط رگرسیون مربوط 
می‌شود و زمانی خودهمبستگی قوی رخ می‌دهد که palis‏ باقیمانده شدیداً باهم در ارتباط باشنده به‌عبارت‌دیگر 
تغییراتشان به‌صورت سیستماتیک رخ دهد. شاخص موران ابزاری در تحلیل خودهمبستگی فضایی است که به تحلیل 
الگوی توزیع عوارض در فضا باملاحظه هم زمان موقعیت مکانی و مقدار خصیصه موردنظر می‌پردازد. نتایج حاصل 
از این تحلیل OLE‏ می‌دهد که آیا عوارض به‌صورت تصادفی, پراکنده یا خوشه‌ای در فضا توزیع شده‌اند. هر چه 
میزان شاخص موران به یک نزدیک‌تر باشد نشان دهنده قوی‌تر بودن الگوی خوشهای است (بلیانی و حکیم 
دوست»۱۳۹۳؛ عسگری»۱۶:۱۳۹۰). 


بر اساس این آزمون شاخحص موران ۰,۹۸۶۱۸۳ است و ازآنجاکه مقدار آن مثبت. بالای صفر و Kay‏ به یک 


است. می‌توان نتیجه گرفت که الگوی E 22 à‏ زمین در دشت مرکزی قائنات دارای خود همبستگی فضایی و 
داری الگوی خوشه‌ای است. با توجه به بالا بودن امتیاز استاندار =Z (Z‏ ۲۸۱,۲۸) و مقدار P-Value‏ صفر است: 


MA‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


می‌توان همبستگی و خوشه‌ای بودن الگوی داده‌ها یا شانعص فرونشست زمین را در منطقه‌ی موردمطالعه تأیید نمود؛ 


که به معنی نامتوازن بودن سطح برداشت آب‌های زیر زمینی و ... است (شکل ۱۱). 


Moran's Index: 0.984189 gl ams cm | mm Critical Valua 

zowcore: 201.204401 mmm rebum) (a scare) 
o.01 یسم‎ - -2.508 

p-value: 0.000000 oos ica asa 1.9546 
9.30 am 1.96 1.65 

- 5 1.62 - 1.65 

0.10 œ 1.65 — 1.96 
0.07 bow 1.94 2 
ooi — - 32^ 


| 


شکل ۱۱- تحلیل الگوی نرخ فرونشست زمین با شاخص خودهمبستگی فضایی در منطقه موردمطالعه 


به‌منظور پیش‌بینی رخداد پدیده فرونشست در آینده از پارامترهای تغییرات سطح آب ارتفاع سطح زمین» شیب و 
ناهمواری‌ها استفاده شد. خلاصه مشخصات آماری پارامترهای موردمطالعه در جدول ۳ نشان oala‏ شده است. نتایج 
جدول نشان می‌دهد که ناهمواری‌های محدوده مطالعاتی دامنه‌ای بین ۱ تا ۳۳ متر داشته و بيشینه شیب منطقه حدود 
YA‏ درجه است. با توجه به برداشت‌های صورت گرفته» تغییرات تراز آب زیرزمینی و میزان افت آن از ۰,۸۲ متر تا 
بیش از ۳۰ متر در منطقه مشاهده می‌شود. میانگین این افت» حدود ۱۶ متر است. بیشترین تراکم مشاهده شده 


فرونشست زمین در منطقه نیز بود. 


جدول ۳. خلاصه مشخصات آماری پارامترهای موردمطالعه 


شاخص کمینه پیشینه میانگین انحراف معیار 
نرخ فرونشست زمین eem - ۹ 0 OF‏ 
سطح برداشت آب زیرزمینی Yow IAS MAYA Wey‏ 
ارتفاع ۱۱۳۳ ۱9۹۷ Yva V‏ ۱۳۹۶۹۹ 
شیب £jev wy 7۸ NE‏ 
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سال یازدهم تحلیل فضایی فرونشست سطح زمین با استفاده از .... ۱۹۹ 


به‌منظور شناخت ارتباط بین این عوامل ماتریس همبستگی محاسبه شد که نتایج آن در جدول ٤‏ نشان داده شده 
است. یافته‌های همبستگی بین شاحص‌ها نشان داد که بین عامل تغییرات تراز سطح آب و نرخ فرونشست زمین» 
همبستگی کمی وجود داشته =p)‏ ۰,۱۳۸-) که این ارتباط غیرمستقیم است؛ یعنی هر چه میزان برداشت آب زیرزمینی 
بیشتر باشد L)‏ سطح تراز اب کمتر شود)» میزان فرونشست نیز بیشتر می‌شود لازم به ذکر است که با توجه به سطح 
معناداری -sig)‏ ۰,۵۸۵) این ارتباط معنادار نیست؛ اما بین نرخ فرونشست زمین و شیب و ارتفاع زمین همبستگی 
معناداری وجود دارد. که با توجه به مقدار آماریه پیرسون» همبستگی بین عوامل» معکوس بوده است؛ به‌عبارت‌دیگر 


هر چه میزان ارتفاع و شیب کمتر شده است. میزان فرونشست زمین بیشتر شده است. 


جدول ۶- ماتریس همبستگی تراکم رخداد فرونشست (متغیر وابسته) با سایر پارامترهای مستقل 


متغیر وابسته نرخ فرونشست زمین 
متغیر مستفل Sig. (2-tailed) Pearson Correlation‏ نتیجه 
سطح برداشت آب زیرزمینی ۳۸ - 2۸۵« عدم ارتباط معنادار 
ارتفاع "£60 - Yl‏ ارتباط معنادار 
شیب eve"‏ - ۷ ارتباط معنادار 
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۵- نتیجه گیری 

اندازه‌گیری میزان و iala‏ جابجایی سطح زمین از طریق اطلاعات ماهواره‌ای روشی نوین در پایش این پدیده است 
که Ub‏ دهه قبل مطرح شده است. تداخل سنجی راداری تا کنون یکی از دقیق‌ترین و کم‌هزینه‌ترین روش‌های 
سنجش‌ازدور برای شناسایی و نمایش جابجایی به وجود آمده در سطح زمین است. همچنین این روش به‌واسطه 
بهره‌گیری از داده‌های ماهواره‌ای و قابلیت تکرارپذیری آن پایش این پدیده در مکان موردنظر با سهولت و کوتاه‌ترین 
زمان امکان اجرا دارد. این روش بررسی میزان جابجایی سطح زمین را در تمام محدودۀ موردبررسی امکان‌پذیر کرده 
است و امکان پایش of‏ را برای تمام نقاط برای دوره‌های زمانی مختلف فراهم نموده است. در این تحقیق از روش 
تداخل سنجی راداری برای تحلیل و اندازه‌گیری میزان جابجایی سطح زمین برای دشت مرکزی شهرستان قائنات 
استفاده شده است. نتایج این بررسی OUS‏ داد که بیشترین فرونشست در بازه زمانی یک‌ساله و فصلی طی سال‌های 
۷ تا ۲۰۲۰ برای منطقه موردمطالعه موردبررسی قرار گرفت. نتایج نشان داد که میزان فرونشست زمین از سال 
۷ تا سال ۲۰۲۰ افزایش داشته است؛ به‌طوری که از ۰,۰۲٩‏ متر در سال ۲۰۱۷ به ۰,۰۵۳۹ متر در سال ۲۰۱۸ 


۲, متر در سال ۹ و در سال ۰ به ۰,۰۷۸ متر VA)‏ سانتی‌متر در سال) می‌رسد. پهنه‌های در معرضص 
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سالانه بیشتر منطبق بر نواحی جنوب شرقی دشت قائنات یعنی نحای پیرامونی شهر اسفدن است؛ اما 


فرو 
نواحی شمال غربی دشت و نواحی پیرامونی شهر قاین نیز با فرونشست همراه است. 

پژوهش‌های پیشین. بر نقش و اهمیت استفاده از تصاویر راداری سنجنده‌های مختلف در شناسایی دقیق این 
پدیده تأکید داشته‌اند. تحقیق حاضر همانند مطالعات امام اوغلو | و همکاران (۲۰۱۸) بای لین" و همکاران MVNA)‏ 
لیو" و همکاران QA)‏ اوسمانو" و همکاران (۲۰۱۳) و کیو" و همکاران YE)‏ از سازگاری بسیار SUL‏ نتایج 
تداخل سنجی راداری با مشاهدات زمینی را OUS‏ می‌دهد و مناسب بودن این روش در اندازه‌گیری تغییر شکل سطح 
زمین را تائید می‌کند. منطقه موردمطالعه در دهه‌های اخیر رشد جمعیت را در کنار توسعه کشاورزی تجربه کرده 
است. بارش کم سال‌های اخیر و خشک‌سالی‌های متوالی. محدودیت منابع آب سطحی, روش‌های نادرست آبیاری 
در بخش کشاورزی در منطقه مطالعاتی» به همراه الگوی کشت نامناسب. موجب حفر بی‌رویه چاه‌های آب در منطقه 
شده است. برداشت بی‌رویه از این چاه‌ها و منابع آب زیرزمینی» موجب عدم تعادل آبخوان منطقه شده و پایین رفتن 


سطح آب زیرزمینی را به دنبال داشته است. 


به‌منظور بررسی ارتباط فرو با برداشت آب زیرزمینی در منطقه, اندازه‌گیری تراز آب. از سطح زمین تا 
سطح آب موجود در چاه انجام شد که هیدروگراف معرف آب زیرزمینی در دشت مرکزی قائنات افت محسوسی در 
ead oh a‏ که SA D‏ ا سنوی ام boss‏ فسال isst OY AA OT dla‏ 
است» یعنی حدوداً سالی ۲٢‏ سانتی‌متر افت سطح آب‌های زیرزمینی را داشته‌ايم. نتایج حاصل از این بخش مطالعه, در 


ت تحت تأثیر 


ارومیه است. تحلیل خودهمبستگی فضایی و شاخص موران نیز خوشه‌ای بودن رخداد فرو 
تغییرات سطح OF‏ در منطقه موردطالعه را تأبید کرد. با توجه به Ub‏ بودن امتیاز استاندار ,۲۸۱,۲۸(2 = (Z‏ و مقدار P-‏ 
Value‏ صفر است. می‌توان همبستگی و خوشه‌ای بودن الگوی داده‌ها یا شاخص فرونشست زمین را در منطقه‌ی 
موردمطالعه تأیید نمود؛ که به معنی نامتوازن بودن سطح برداشت آب‌های زیر زمینی و ... است. نتایج همبستگی بین 
شاخص‌ها نشان داد که بین عامل تغییرات تراز سطح آب و نرخ فرونشست زمین» همبستگی کمی وجود داشته 
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سال یازدهم تحلیل فضایی فرونشست سطح زمین با استفاده nel‏ ۱۳ 


(p = -۰,۱۳۸(‏ که این ارتباط غیرمستقیم است؛ یعنی هر چه میزان برداشت آب زیرزمینی بیشتر باشد L)‏ سطح تراز 


آب کمتر شود) میزان فرونشست نیز بیشتر می‌شود لازم به ذکر است که با توجه به سطح معناداری =sig)‏ ۰,۵۸۵) 
این ارتباط معنادار نیست که با یافته‌های حاصل از به‌کارگیری این روش با دیدگاه افت سطح آب زیرزمینی که توسط 
پژوهشگران مختلف مانند گالووی" و همکاران (۲۰۱۱) متاگ" و همکاران CHA)‏ فین " و همکاران (۲۰۰7) همسو 
نیست؛ اما نتایج به دست آمده با تحقیقاتی که توسط آروین و همکاران (YA)‏ حاجب و همکاران (PAA)‏ دهقانی 
“lol YAE)‏ و همکاران X)‏ شوجون* و همکاران QVO‏ توماس" و همکاران (۲۰۱۰)» Sigg‏ " و همکاران 
(T+ +4)‏ انجام شد مطابقت دارد. همچنین بر اساس نتایج پژوهش حاضر بین نرخ فرونشست زمین و شیب و ارتفاع 
زمین همبستگی معناداری وجود دارد. که با توجه به مقدار آماریه پیرسون» همبستگی بین عوامل. معکوس بوده است؛ 
به‌عبارت‌دیگر هر چه میزان ارتفاع و شیب کمتر شده است. میزان فرونشست زمین بیشتر شده است. 

لازم به یاداوری است که کاهش بارندگی» بروز حشکسالی‌های پیوسته در چند سال اخیر و برداشت‌های بی‌رویه 
و غیرمجاز» در کاهش حجم ذخایر سفره و بروز روند منفی عمق دسترسی به آب زیرزمینی آبخوان دشت قائنات 
بی‌تأثیر نبوده است. همچنین در دهه‌های گذشته روند منفی برای تراز آب زیرزمینی در نقاط مختلف جهان و ایران 
گزارش شده است که با توجه به ld‏ اقلیمی خشک و نیمه‌عشک کشور که Ol ul‏ را کشوری دارای تنش آبی زیاد 
معرفی کرده است. 

درنتیجه» این محدوده نه‌تنها ازنظر کمی در وضعیت نامناسبی قرار دارد. ازنظر کیفی هم در معرض خطر خواهد 
بود. به همین دلیل» برای کنترل این پدیده و کاهش اثرات منفی OT‏ به مدیریت درست و برنامه‌ریزی مناسب‌تر نیاز 
خواهد بود؛ و اگر برداشت بی‌رویه آب به همین ترتیب ادامه یابد بدون شک در آینده عواقب نامطلوب در مورد منابع 
آب زیرزمینی دشت قائنات و به‌تبع آن بحران‌های اجتماعی اقتصادی و سیاسی در منطقه اتفاق خواهد افتاد. تغذیه 
مصنوعی سفره آب زیرزمینی دشت مرکزی قائنات و تجدیدنظر در میزان بهره‌برداری از آبخوان مربوطه. جلوگیری از 
حفر چاه‌های غیرمجاز, افزایش راندمان آبیاری» جلوگیری از کشت محصولات نیازمند به آب فراوان» پوشش انهار و 
صرفه‌جویی در مصرف آب و آموزش کشاورزان می‌تواند تا حدودی از افت بی‌رویه سطح آب این دشت جلوگیری 


شود. 
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